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The spins of a number of supermassive and stellar-mass black holes have been measured based 
on detections of thermal continuum emission and relativistically broadened iron lines in their x-ray 
spectra. Likewise, quasiperiodic variability has been observed in several sources. Such measurements 
commonly make the assumption that black holes are described by the Kerr metric, which according 
to the no-hair theorem characterizes black holes uniquely in terms of their masses and spins. This 
fundamental property of black holes can be tested observationally by measuring potential deviations 
from the Kerr metric introduced by a parametrically deformed Kerr-like spacetime. Thermal spectra, 
iron lines, and variability have already been studied extensively in several such metrics, which 
usually depend on only one particular type of deviation or contain unphysical regions outside of 
the compact object. In this paper, I study these x-ray probes in the background of a new Kerr-like 
metric which depends on four independent deviation functions and is free of pathological regions 
outside of the event horizon. I show that the observed signals depend significantly on primarily two 
types of deviations and that the strong correlation between the spin and the deviation parameters 
found previously in other Kerr-like metrics is partially broken for rapidly spinning black holes. This 
suggests that high-spin sources are the best candidates for tests of the no-hair theorem with x-rays 
and I obtain first constraints on such deviations from the stellar-mass black hole Cygnus X-1. 

PACS numbers: 04.50.Kd,04.70.-s 


I. INTRODUCTION 

Black holes are encountered ubiouitously in the centers 
of galaxies and in x-ray binaries [l|, l2| . In general relativ¬ 
ity, stationary black holes in vacuum are fully character¬ 
ized by the no-hair theorem which states that they only 
depend on two parameters, their masses M and spins J, 
and are described by the Kerr metric (see, e.g.. Ref. i). 
In practice, such black holes will not be fully stationary 
or exist in perfect vacuum due to the presence of, e.g., 
surrounding matter or other stars, but it is often reason¬ 
able to assume that these effects are sufficiently small to 
be negligible. While observational evidence exists for the 
presence of event horizons (e.g., i), a final proof of the 
Kerr nature of these compact objects is still lacking. 

General relativity has been well confirmed in the 
regime of weak spacetime curvatures @ but still remains 
practically untested in the strong-held regime that is 
found around compact objects (see Ref. @). Indeed, 
black holes need not necessarily be described by the Kerr 
metric at all (e.g., 0)- Since it is unclear at present 
whether general relativity is modihed in the strong-held 
regime and, if so, in what manner, an efficient approach 
is to test the no-hair theorem in a model-independent 
manner that is based on a parametrically deformed Kerr- 
like metric which encompasses many different theories of 
gravity at once. Such metrics generally do not derive 
from the action of any particular gravity theory, and the 
underlying theory is usually unknown. However, it is 
assumed that particles follow the geodesics of the space- 
time. 

Several such metrics have been constructed thus far 


(e.g., [Ml), which depend on one or more free param¬ 
eters that measure potential deviations from the Kerr 
metric and which include the Kerr metric as the special 
case when all deviations vanish. Observations can then 
be used to measure these deviations, should they exist, 
and, thereby, infer properties of the underlying but gen¬ 
erally unknown theory of gravity. If no deviations are de¬ 
tected, the compact object is verified to be a Kerr black 
hole. If, on the other hand, nonzero deviations are mea¬ 
sured, there are two possible interpretations. If general 
relativity still holds, the object is not a black hole but, 
instead, another stable stellar configuration or, perhaps, 
an exotic object M- Otherwise, the no-hair theorem 
would be falsified. 

Tests of the no-hair theorem have been suggested based 
on either the gravitational-wave signals of extreme mass- 
ratio inspirals (see Ref. [Tsj l for a review) or electromag¬ 
netic observations of pulsar black hole binaries [l^ , stars 
on orbits around Sgr A* [13 , or of the radiation from 
the accretion flows of black holes (e.g., M)- The emis¬ 
sion from black hole accretion disks provides a number of 
different observables that could be used for such tests, in¬ 
cluding the direct imaging of the shadows of supermassive 
black holes with very-long baseline interferometric arrays 
(e.g.. HOl - l^ b relativistically broadened iron lines (e.g., 
[23l - l25j b thermal disk spectra (e.g., quasiperi- 

odic variability , and x-ray polarization . 

Previous model-independent studies in the electromag¬ 
netic spectrum were largely based on either a quasi-Kerr 
metric or a metric proposed by Ref. [13 • These met¬ 
rics depend on only one type of potential deviations from 
the Kerr metric, although there are at least four inde- 
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pendent ways to parametrize them (cf. El)- Both of 
these metrics generally harbor naked singularities which 
makes them less favorable candidates as an alternative 
to black holes. In addition, the quasi-Kerr metric, just 
as many other Kerr-like metrics, contains pathological 
regions near the central object where other singularities 
occur or causality is violated. Therefore, an artificial 
cutoff needs to be imposed in both metrics which shields 
the singularities and any unphysical region from outside 
observers. In this form, then, these metrics are valid 
frameworks for tests of the no-hair theorem (see Ref. [32| 
for a full discussion). Note, however, that the metric of 
Ref. [l^ uses the Newman-Janis algorithm [s^ to gen¬ 
erate rotating solutions from static seeds and it is not 
guaranteed that this procedure can be applied consis¬ 
tently to general metrics which are not solutions of the 
Einstein field equations. 

Under the assumption that black holes are described 
by the Kerr metric, thermal disk spectra and relativisti- 
cally broadened iron lines have already been widely used 
to measure the spins of stellar-mass black holes and ac¬ 
tive galactic nuclei (AGNs); see Refs, (s^, [s^ for recent 
reviews. The latter method assumes irradiation off the 
accretion disks of black holes in the form of fluorescent 
iron lines which appear significantly broadened due to the 
relativistic effects of light bending, Doppler boosting and 
beaming, and the gravitational redshift [^[s^. The for¬ 
mer method is based on the Novikov-Thorne model [3^ . 
a relativistic thin-disk accretion flow model, which de¬ 
scribes the emission of thermal radiation. Several authors 
have computed the observed spectra with the inclusion 
of the relativistic effects |39l44ll| . The current state-of- 
the-art fitting routine combines such spectra with a suit¬ 
able model of the disk atmosphere to account for spectral 
hardening [d^. Both of these methods aim to measure 
the location of the innermost stable circular orbit (ISCO) 
which depends directly on the spin of the black hole. 

Quasiperiodic oscillations (QPOs) have been observed 
in both stellar-mass black holes and AGNs [i, 113 and 
provide a potential third means to measure their spins. 
QPOs are transient and occur mostly in nonthermal ac¬ 
cretion disk states and during state transitions. They can 
be divided into two general classes: high-frequency QPOs 
(roughly 40-450 Hz) and low-frequency QPOs (roughly 
0.1-30 Hz) and their observed frequencies sometimes fall 
into ratios of small integers; see Ref. 0 for a detailed dis¬ 
cussion. A number of different QPO models have been 
designed to date (e.g., [43 - l48l| l. but at present the exact 
origin of QPOs is still unknown. 

In this paper, I analyze thermal disk spectra, relativis- 
tically broadened iron lines, and quasiperiodic variability 
in the context of a new Kerr-like metric [d^. This met¬ 
ric describes a black hole and is free of singularities or 
other adverse regions outside of the event horizon. Like¬ 
wise, this metric depends on four independent deviation 
functions which parametrize possible departures from the 
Kerr metric. I show that these observables depend pri¬ 
marily on two deviation types which can lead to signif¬ 


icantly modified signals that should be resolvable with 
either current or near-future x-ray missions. 

While the spin can be inferred uniquely from the loca¬ 
tion of the ISGO for Kerr black holes, in Kerr-like metrics 
the ISGO generally depends on the spin and the deviation 
parameters (e.g., 0). References [1^ [H, showed 
that iron line profiles are practically indistinguishable for 
values of the spin and the deviation parameter that cor¬ 
respond to the same ISGO in either the quasi-Kerr met¬ 
ric 0 or the metric of Ref. 0|. Reference found 
a similar result for thermal disk spectra obtained in the 
latter metric. Thus, there is a strong correlation between 
the spin and the respective deviation parameters which 
hampers the ability of these observables (as well as of 
quasiperiodic variability) to clearly measure such devia¬ 
tions if they exist. 

I demonstrate that both the thermal spectra and iron 
line profiles are also very similar in the new Kerr-like met¬ 
ric for values of the spin and the deviation parameters 
with the same ISGO if the ISGO is located at compara¬ 
tively large radii corresponding to Kerr black holes with 
small to intermediate spins. However, for ISCO radii that 
are sufficiently close to the event horizon, i.e., for rapidly 
spinning black holes, this strong correlation among the 
spin and the deviation parameters is partially broken and 
the thermal spectra and line profiles become distinguish¬ 
able. I also show that a similar property holds for the 
QPOs predicted by the diskoseismology (see Ref. 0 for 
a review) and resonance [U 0 models. 

This suggests that black holes with high spins are the 
best candidates for tests of the no-hair theorem with x- 
ray observables. Based on the recent near-maximal spin 
measurement of the Galactic black hole Cygnus X-1 from 
its thermal spectrum by Gou et al. [0 . 1 obtain first con¬ 
straints on two deviation parameters of the new metric. 

In Sec. im I summarize some of the properties of the 
Kerr and Kerr-like metrics. In Sec. imi I describe the 
ray tracing algorithm that I use to simulate the observed 
spectra and line profiles. In Secs. Hvl El and ED I de¬ 
scribe the effects of the deviation parameters on the ther¬ 
mal spectra, iron lines, and QPOs, respectively. I discuss 
my conclusions in Sec. lVIll From here on, I use geometric 
units and set G = c = 1, where G and c are the gravita¬ 
tional constant and the speed of light, respectively, unless 
otherwise stated. In these units, length (and time) scales 
are expressed in terms of the mass of the black hole, 
where M = GM/c^ = Vg, the gravitational radius. 


II. KERR AND KERR-LIKE BLACK HOLE 
METRICS 

In this section, I describe some of the properties of 
black holes that are described by the two metrics I con¬ 
sider in this paper. Expressed in Boyer-Lindquist coor- 
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dinates, these are the Kerr nietric 
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In these expressions, a = J/M is the spin parameter. 

The latter metric contains the four free functions /(r), 
Ai(r), A 2 {r), and ^ 5 (r) that depend on four sets of pa¬ 
rameters which measure potential deviations from the 
Kerr metric. In the case when all deviation parameters 
vanish, i.e., when /(r) = 0, Ai{r) = A 2 {r) = A^{r) = 1, 
this metric reduces to the Kerr metric in Eq. CD- 

Just as the Kerr metric, the metric in Eq. © is sta¬ 
tionary, axisymmetric, and asymptotically flat. Unlike 
the Kerr metric, which is a vacuum solution of the Ein¬ 
stein equations, the metric in Eq. is not a vacuum 
solution in general relativity, because it describes black 
holes in modified theories of gravity that differ from Kerr 


black holes. Thanks to its structure, the Kerr-like metric 
in Eq. © is phenomenological; i.e., it encompasses large 
classes of modified gravity theories and can be mapped 
to the known black hole solutions of particular theo¬ 
ries for certain choices of the deviation functions includ¬ 
ing Chern-Simons gravity and Einstein-Dilaton-Gauss- 
Bonnet gravity (see the appendix of Ref. [i^ for the re¬ 
spective mappings). Using this metric as a framework, 
one can then test observationally whether astrophysical 
black holes are indeed described by the Kerr metric. 

In general relativity, stationary, axisymmetric, and 
asymptotically flat metrics that admit the existence of in- 
tegrable two-dimensional hypersurfaces generally depend 
on only four functions. As a consequence of Frobenius’s 
theorem, such hypersurfaces are automatically guaran¬ 
teed to exist if the spacetime is also vacuum. Such met¬ 
rics can be written in the form of the Papapetrou line 
element where the metric is expressed with respect to the 
Weyl-Papapetrou coordinates and has only three metric 
functions. Out of these functions only two are indepen¬ 
dent, while the third one can be derived from the other 
two. 

This happens for three reasons. First, the symmetries 
imposed, the assumption of asymptotic flatness and the 
vanishing of the Ricci tensor allow for the spacetime to 
have integrable two-dimensional hypersurfaces that are 
orthogonal to the two Killing fields and on which one 
can define coordinates that can be carried along integral 
curves of these Killing fields to the rest of the spacetime. 
Thus the metric can be written in a 2 x 2 block form. 
If one chooses as one of the coordinates the determinant 
of the {t, (f) part of the metric, then the metric can be 
written in a form that has only four independent func¬ 
tions. Second, the field equations imply, by the vanishing 
of the Ricci tensor (i.e., the vacuum assumption), that 
the coordinate p which is defined by the determinant of 
the (f, (j)) part of the metric is a harmonic function and 
thus one can define the second coordinate on the two- 
dimensional hypersurfaces as the harmonic conjugate of 
p and absorb one of the functions in the process, reducing 
the independent functions to three. Finally the vacuum 
field equations imply that the third of the functions is 
related to the other two and can be determined up to the 
addition of a constant [^ . 

In alternative theories of gravity, however, the vac¬ 
uum assumption does not necessarily imply the existence 
of two-dimensional integrable hypersurfaces. Therefore, 
one would expect that metrics which describe black holes 
in alternative theories of gravity depend on at least four 
independent functions which motivates the introduction 
of four deviation functions in Eq. ( 151 ) . The quasi-Kerr 
metric M and the metric of Ref. |l2 | depend on only 
one such function making them less general in this sense. 
Nonetheless, these metrics can likewise be used for ob¬ 
servational tests of the nature of black holes. Recently, 
Cardoso et al. generalized the metric of Ref. [l2| to 
include two deviation types. 

The deviation functions in Eqs. ©-© are written as 
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power series in M/r. The lowest-order coefficients of 
these series vanish so that the deviations from the Kerr 
metric are consistent with all current weak-field tests of 
general relativity (see Ref. i). In this paper, I will fo¬ 
cus on black holes that are described by the lowest-order 
metric, which only depend on the mass M, spin a, and 
the deviation parameters £ 3 , 013 , 022 , and a 52 . 

Certain lower limits on the deviation parameters £ 3 , 
ai 3 , 0 ^ 22 , and 052 exist which are determined by the prop¬ 
erties of the event horizon. In order for an event horizon 
to be present these parameters must obey the relations 


where 


£3 > B'i, ai3 > Rs, 
022 > I?2, 052 > R 2 , 


B2 = 

B3 = 


{M + 
{M + 
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For values of the deviation parameters smaller than these 
limits the metric harbors a naked singularity instead of 
a black hole [i^ . The event horizon itself is independent 
of these parameters and coincides with the event horizon 
of a Kerr black hole. 


r+ = M+y/M^ - a?-. (11) 

As discussed in Refs. [ll|,|i^ the deviation parameters af¬ 
fect the locations of the circular photon orbit and ISCO, 
as well as the energy, axial angular momentum, and dy¬ 
namical frequencies of particles on circular equatorial or¬ 
bits. For nonzero values of the deviation parameters the 
location of the ISCO is shifted and the orbital velocity of 
the disk particles either increases or decreases [i^ (cf. 
Refs. [H, [ 1 ^). In addition, photons which are emit¬ 
ted by an accretion disk surrounding the black hole and 
which are detected by a distant observer also experience a 
modified gravitational redshift, relativistic boostin g an d 
beaming, as well as light bending (cf. Refs. [I 3 , 1^1. 
These properties can have a si gnifi cant effect on the ob¬ 
served spectra and QPOs (see PHI). 

Both the Kerr and the Kerr-like metrics admit three 
constants of motion thanks to their symmetries . 

This means that for a particle with rest mass ^ and 4- 
momentum p°‘ on a geodesic orbit its energy E, angular 
momentum about the rotation axis of the black hole, 
and the so-called Carter constant are conserved. The 
Carter constant has the same form in both metrics, given 
by the expression 

Q = p0 — {Lz—aE)'^+pL^a^ cos^ 9-\ - \—{L^—aE sin^ 0)^. 

sin 6 

( 12 ) 

While in general relativity the no-hair theorem estab¬ 
lishes that the Kerr metric is the unique black hole metric 
(assuming the black hole is electrically neutral), in alter¬ 
native theories of gravity a Kerr-like metric generally pos¬ 
sesses only two constants of motion which correspond to 


its stationarity and axisymmetry. Therefore, the choice 
of the metric in Eq. ([3]) necessarily restricts potential de¬ 
viations from the Kerr metric to belong to the class of 
Kerr-like metrics with three constants of motion and it 
would be desirable to employ an even more general Kerr- 
like metric which contains all Kerr-like metrics with at 
least two constants of motion. At present, however, no 
such metric is known. 

On the other hand, the existence of twodimentional in- 
tegrable hypersurfaces of the metric in Eq. o has two 
distinct advantages. Eirst, the known black hole solu¬ 
tions of alternative theories of gravity (all of which can 
be mapped to this metric [i^) possess three constants of 
motion at least in the case of slowly rotating black holes 
[ 5 ^ . Second, the existence of three constants of motion 
is of great practical use and simplifies many analytic and 
numerical computations. I will make use of this property 
in the following. 


III. RAY TRACING ALGORITHM 


The existence of three constants of motion can sim¬ 
plify the ray tracing computations required for the study 
of the emission from the accretion disks of black holes. 
Thanks to these constants, all four photon equations of 
motion separate and can be written in first-order form. 
Eor the Kerr-like metric in Eq. ©, these are given by the 
equations [i^ 


= (13) 


, dr 




dX sin^ 0 
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where 


P(r) = (r^-b a^)Ai(r) - aA 2 (r)^, (17) 

R(r) = P(r)2-A[(C-a)2+ry], (18) 

0(0) = ry-b cos^ 0 — cot^ 0. (19) 


In these expressions, A is an affine parameter, and I used 
the standard abbreviations 


^ E ’ 

(20) 

Q 
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Equivalently, these equations can be written in an in¬ 
tegral form, which contains only elliptic integrals in the 
case of the Kerr metric. These can then be calculated 
elegantly and fast [13, [13|- Foi' the Kerr-like metric in 
Eq. (|3]), however, the equations of motion are signif¬ 
icantly more complicated and cannot be expressed in 
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terms of simple elliptic integrals. Therefore, the inte¬ 
gral form of these equations does not lead to an increase 
in computation speed and I integrate the equations of 
motion instead. 

However, a direct integration of Eqs. (IT^-(IT6l) suf¬ 
fers from an unacceptable loss of precision around the 
radial and polar turning points where the radial or po¬ 
lar momenta vanish, respectively. Therefore, I follow the 
scheme used by Ref. who solved the first-order equa¬ 
tions for the coordinates t and 4> and the second-order 
geodesic equations for r and 9. 

As in Ref. , I consider an accretion disk that lies in 
the equatorial plane of the black hole with a radial extent 
from the ISCO to some outer radius rout and that is 
viewed by an observer at an inclination angle * at a large 
distance ro from the black hole. This observer defines 
an image plane, where photons emitted from the disk are 
observed. This image plane can be set up in the local rest 
frame of the observer and has the Cartesian coordinate 
axes (^.[ 5 ^ 


X =-ro - y - w, 

(1 + 

/ ±\/0(o 

V =ro - 7 - 

v^< (1 + Sff) 


( 22 ) 

(23) 


In the limit rp —t 00 , these expressions reduce to the 
simple form (cf. [Hi) 



sinz ’ 


y' = 


(24) 

(25) 


Note that the y' axis is chosen in a manner such that it 
lies in a plane with the rotation axis of the black hole 
and that the origin of the image plane is centered on the 
black hole. 

Since the quantities ^ and 77 are constant along the 
photon trajectories, I can calculate them at any given 
point. For each photon corresponding to a point {x',y') 
in the image plane, I invert Eqs. (1^^ and (1^51) to obtain 
the corresponding constants of motion ^ and y and calcu¬ 
late the initial location and 4-momentum of the photon 
expressed in the coordinate frame centered on the black 
hole (cf. [13 )• I then integrate the equations of motion 
(inii-dini) backwards in time until the photon reaches the 
equatorial plane, where I specify the appropriate bound¬ 
ary conditions for the emissivity depending on the appli¬ 
cation as I describe in Secs. m and IS 

The integrability of the equations of motion can also 
be used to calculate analytically the observed size of the 
accretion disk in the image plane. The equations of mo¬ 
tion, then, only have to be integrated for image points 
that lie inside of this region, because only photons cor¬ 
responding to these points reach the equatorial plane in 
the radial interval between the ISCO and the outer disk 
radius. Since the functions R{r) and 0(0) depend of the 
constants of motion ^ and y, this set of image points can 


be obtained from the requirement that these functions 
be non-negative. These points form a disk in the image 
plane, which for large disks extends to a radius that is 
only slightly larger than rout (cf. [H3)- Thanks to the 
gravitational attraction of the black hole, it is likewise 
sufficient to integrate the equations of motion only for 
photons with image points at a radius that is larger than 
the ISCO. 

IV. THERMAL DISK SPECTRA 

In this section, I summarize the relevant properties of 
the thin-disk accretion flow model that I use and calcu¬ 
late the observed thermal disk flux density. I proceed to 
analyze the effects of the deviation parameters on both 
the emitted and the observed flux densities. 


A. Thin accretion disk model 


For the accretion disk of the black hole, I assume the 
standard relativistic thin-disk model of Ref. [s^ . In this 
model, the accretion disk lies in the equatorial plane of 
the black hole and is geometrically thin and optically 
thick. It extends from the ISCO to some outer radius 
Tout) and the disk particles move on nearly circular equa¬ 
torial orbits as they are accreted by the black hole. In 
my analysis, I likewise assume that there is no torque at 
the inner boundary of the disk which is located at the 
ISCO. 

Using the conservation of rest mass, energy, and an¬ 
gular momentum of the disk particles. Ref. derived 
the time-averaged thermal energy flux emitted from an 
accretion disk. The result of Ref. is valid for all 
black holes that are stationary, axisymmetric, asymptot¬ 
ically flat, and reflection-symmetric across the equatorial 
plane. Therefore, their model applies not only to the 
Kerr metric, but also to the Kerr-like metric in Eq. (1^. 

The emitted energy flux is given by the expression [60| 


where M is the mass accretion rate and 
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is a dimensionless function of radius. In the last expres¬ 
sion, g is the determinant of the (t, r, (j)) part of the met¬ 
ric. For the Kerr metric, -^/—g = r, while for the metric 
in Eq. dS]) 
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FIG. 1. Radial profiles of the luminosity of an accretion disk 
(in units of c®/ 47 rG^) around a black hole with spin a = 0.9M 
for several values of the deviation parameters ais (top) and 
ce 22 (bottom). In each panel, only one deviation parameter is 
varied, while the other one is set to zero. The luminosity pro¬ 
file is strongly peaked at a radius near the ISCO and decreases 
rapidly at larger and smaller disk radii. For decreasing values 
of the parameters ais as well as for increasing values of the 
parameter 0 : 22 , the peak of the luminosity profile increases 
and is located closer to the ISCO radius. 


From the energy flux in Eq. (PSI) . I obtain the disk 
luminosity 


L = 2 f 2'KrF{r)dr. (29) 

In Fig. [1] I plot the radial profile of the disk luminosity 
for a black hole with spin a = 0.9M for different values 
of the deviation parameters ais and a 22 - The luminosity 
profile has a sharp peak at a radius near the ISCO and 
falls off quickly at larger and smaller disk radii. At the 


ISCO, the emitted flux is zero. As either the value of the 
parameter aia decreases or the value of the parameter 
a 22 increases, the peak of the luminosity profile increases 
and is located at a radius that lies closer to the ISCO. The 
deviation parameters 63 and q ;52 have only a minor effect 
on the luminosity profile which primarily either increases 
or decreases its maximum slightly. 

Since the ISCO in the Kerr-like metric in Eq. (jS]) de¬ 
pends on the spin and on the deviation parameters ais 
and a 22 (as well as to a lesser extent on the parameter 63 ), 
it is important to analyze the dependence of the luminos¬ 
ity profile on different combinations of these parameters 
which correspond to the same ISCO. In Fig. [2 I plot 
several cases of luminosity profiles for Kerr black holes 
with low, intermediate, and high spins together with the 
luminosity profiles for Kerr-like black holes with different 
spins and nonzero values of either the parameter ais or 
a 22 which match the ISCO of the respective Kerr black 
hole. 

The effect of both deviation parameters is very similar. 
For small to intermediate values of the spin, the lumi¬ 
nosity profile changes only slightly even for large values 
of the deviation parameters, because the ISCO and the 
peak of the luminosity profile are located at relatively 
large radii. For high spin values, however, the luminosity 
profiles are altered significantly. 

Since the disk is in thermal equilibrium, the disk flux is 
related to the disk temperature by the standard Stefan- 
Boltzmann law, 


F(r) = ar4,(r), 
with the specific intensity 

2hiy^ 1 T 

^eiye) — 


/< 


( - 


(30) 


(31) 




In these expressions, Teff is the effective temperature, 
Tcoi(r) = /coiT'eff(r) is the color temperature with spec¬ 
tral hardening factor /coi, Ve is the frequency of emission, 
and h and ks are the Planck and Boltzmann constants, 
respectively. Photons emitted from the disk and observed 
at a large distance from the black hole experience a red- 
shift 

_ Apbs _ g;jt/,obsfeobs'^obs /oo'i 

Ee gtj.u,ekeK 

where the subscripts “obs” and “e” refer to the image 
plane and the accretion disk, respectively, and where 
is the Keplerian frequency of the disk particles calculated 
in Ref. 0. 

The factor T in Eq. (1511) accounts for the emission type 
which I assume to be either isotropic or limb darkened. 
I use a standard limb darkening law (see Refs. 00 ) 
which modifies the specific intensity of the disk by the 
factor 


T(re) 



isotropic radiation 
I cos C limb darkened radiation. 


(33) 
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FIG. 2 . Radial profiles of the emitted disk luminosity (in units of c®/47rG^) of black holes with different values of the spin and 
the deviation parameters 0:13 and 0:22 such that the parameter combinations in each panel correspond to the same ISCO. The 
profiles are very similar for sufficiently large ISCO radii corresponding to Kerr black holes with small to intermediate spins. At 
high spin values, however, the emitted fluxes differ significantly. 


where re is the radius of emission. In this expression, 
C is the emission angle relative to the disk normal as 
measured in a local inertial frame. This angle is given by 
the equation 


cosC = 


gV^ 



(34) 


where 77 is given by Eq. (ED. I derive this expression for 
the emission angle in the Appendix. 


B. Observed spectra 


In order to calculate the observed disk flux density, I 
employ the ray tracing algorithm described in the previ¬ 
ous section. Following Ref. [4l|, I calculate the observed 
photon number flux density = FE„t,jEohs, which 

is given by the integral 


^Soba = 


Fobs 


= Nr 


[ 7obs(^obs)dfIobs — jb [ g -^e(^e)dfIobs 
J n/obs J 

f 1 Tr'dr'dip' 


f Fobs 


M 2 


exp 


Ni 

gFU< 


( £^nhs \ 

V keV / 


- 1 


,(35) 


where EJ 


iVo = 0.07205/-f 





N, ^ 0.1331/-J 


M 


-1/4 


10 ^® g s 


-1 


/ M 


(37) 


and where the Lorentz invariant J/i/® was used in 
Eq. ((551) to relate the observed specific intensity to the 
emitted specific intensity. I choose polar coordinates 
in the image plane and an outer disk radius 


?'out = 10®M. I evaluate the integral in Eq. (I35|) numer¬ 
ically using a large number of rays with logarithmically 
spaced radii and equally spaced polar angles in the image 
plane. 

For a Kerr black hole, the observed spectrum of the 
number flux density was discussed in great detail by, 
e.g., Ref. [dlj. For higher values of the spin a of the 
black hole, the spectrum becomes harder, because the 
inner edge of the accretion disk extends to smaller radii, 
which causes the radiative efficiency and the disk temper¬ 
ature to increase. For larger values of the disk inclination 
i, the observed flux density at lower energies decreases, 
because it originates at larger disk radii, where the ob¬ 
served flux density is oc cos I, while the observed flux 
density increases at higher energies, because it originates 
at smaller radii, where the relativistic effects of boosting, 
beaming, and light bending become important. As ex¬ 
pected, when the accretion rate increases, the observed 
flux density likewise increases, and the spectrum becomes 
harder for hig her values of the spectral hardening factor. 
Reference also studied the effects of a nonzero torque 
at the inner edge of the disk and of returning radiation 
(cf. [13). They showed that these effects can be compen¬ 
sated by an adjustment of the mass accretion rate and 
the spectral hardening factor in a disk model without 
returning radiation and a torque at the ISCO. The ef¬ 
fect of limb darkening is particularly strong at high disk 
inclinations and leads to a lower observed flux density 
relative to the observed flux density when the emission 
is isotropic. 


I verified these spectral properties for Kerr black holes 
using the algorithm described in Sec. IIIII In Fig. |51 I 
compare one spectrum generated by this algorithm with 
a spectrum from Ref. [4l| with the same parameters and 
find excellent agreement. Note, however, that my algo¬ 
rithm does not incorporate the effects of returning radia¬ 
tion or a nonzero torque at the ISCO. While these effects 
can be absorbed into modified values of the mass accre¬ 
tion rate and the spectral hardenin g fa ctor in the case of 
a Kerr black hole as shown by Ref. ^ , it is not obvious 



















FIG. 3. Observed number flux density from a geometri¬ 
cally thin Novikov-Thorne type accretion disk around a Kerr 
black hole with spin a — 0.999M, inclination i = 30°, mass 
M = IOM 0 , distance D = 10 kpc, mass accretion rate 
M = 10^® g s“^, and spectral hardening factor fcoi = 1- 
The red triangles denote the number flux density reported by 
Ref. [ 4 ^, while the solid blue line is the number flux density 
calculated with the algorithm described here. 


that this adjustment works equally well in the extended 
theory where the black hole is described by the metric in 
Eq. ([3]), because the deviation parameters may have dif¬ 
ferent effects on the light bending and the orbital velocity 
of the disk plasma at the ISCO. Nonetheless, one would 
expect that such an adjustment of the mass accretion 
rate and the spectral hardening factor can compensate 
for these effects at least approximately, especially when 
the deviations from the Kerr metric are small. 

The effect of the deviations from the Kerr metric on 
the observed spectra depend on the particular parameter. 
The spectra are strongly affected by both the parameters 
013 and 022 in a manner that is similar to the effect of 
changing the spin of the black hole. In Fig. |4l I plot 
the observed spectrum of the number flux density for a 
black hole with mass M = IOMq and spin a = 0.8M 
for different values of the deviation parameters 013 and 
022 . When one deviation parameter is varied, the other 
one is set to zero. The other system parameters are the 
inclination i = 30°, the distance D = 10 kpc, the mass 
accretion rate M = 10^® g s“^, and the spectra hard¬ 
ening factor /col = 1.7, and I assume an isotropic disk 
emission. For decreasing values of the parameter 013 and 
increasing values of the parameter 022 , the number flux 
density extends to higher energies, while the low-energy 
part of these spectra is practically unaffected. Since the 
deviation parameters predominantly affect the inner part 
of the disk, this is as expected. Changing the disk incli¬ 
nation if one or more of the deviation parameters are 



Energy (keV) 


FIG. 4. Observed thermal spectra from a geometrically thin 
accretion disk around black holes with spin a = 0.8M, incli¬ 
nation i = 30°, mass M = IOMq, distance D = 10 kpc, mass 
accretion rate M — 10^^ g and spectral hardening factor 
/col = 1.7 for different values of the deviation parameters 013 
(top) and Q 22 (bottom). The spectra become harder for de¬ 
creasing values of the parameter 013 and increasing values of 
the parameter 0 : 22 . For all spectra, an isotropic disk emission 
is assumed. 


nonzero has qualitatively the same effect as in the case 
of a Kerr black hole. 

Since the emitted flux depends only marginally on the 
parameters €3 and 052 , the effect of these parameters on 
the observed spectra is very small, at least in the case of 
isotropic emission. Both parameters primarily affect the 
peak flux density of the emitted radiation and the over¬ 
all normalization of the observed flux density through 
the determinant of the metric in Eq. (ESI), while nonzero 
values of the parameter €3 also cause a slight shift of the 
location of the ISCO (see Fig. 6 in Ref. j^). Since the 
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FIG. 5. Observed spectra of the number flux density from a geometrically thin accretion disk around black holes with different 
values of the spin and deviation parameters aia and 0:22 such that the ISCO coincides for the parameter combinations shown 
in each panel. The other parameters are the inclination i = 30°, the mass M = IOMq, the distance D = 10 kpc, the mass 
accretion rate M = 10^® g s~^, and the spectral hardening factor fcoi ~ 1.7. The spectra differ significantly at the high-energy 
end, but some of the spectra of the Kerr-like black holes can still be confused with the spectra of a Kerr black hole with a 
different ISCO and spin. In the left panel, the spectrum of the Kerr-like black hole is very similar to the spectrum of a Kerr 
black hole with a spin a ~ 0.96M (not shown). In the center panel, the spectrum of the Kerr-like black hole can barely be 
mimicked by the spectrum of a Kerr black hole with a spin a = 0.999M (red dotted curve). In the right panel, however, the 
spectrum of the Kerr-like black hole cannot originate from a Kerr black hole even with a spin a — 0.999M (red dotted curve), 
because it does not extend to such high energies. 


normalization of the flux density is used in practice to 
infer the mass accretion rate or, equivalently, the disk lu¬ 
minosity in units of the Eddington luminosity (see, e.g.. 
Ref. 0), the parameters £3 and 052 cannot be obtained 
from the flux density normalization. Conversely, since 
their effect on the normalization is negligible, the mass 
accretion rate can be Inferred robustly even If these two 
parameters are nonzero. 

In the case of limb darkened emission, the deviation pa¬ 
rameters £3, ai3, and a 22 have an effect on the observed 
spectra that is comparable In magnitude and analogous 
to the effect of limb darkening for Kerr black holes (cf. 
(4ll|). As before, the parameter 052 has only a very mi¬ 
nor effect on the observed spectra. Consequently, the 
observed flux density depends strongly on the deviation 
parameters ai 3 and 022 in the case of isotropic emis¬ 
sion, while the observed flux density depends strongly on 
the deviation parameters £3, ais, and 022 in the case of 
limb darkened emission. In all cases, these parameters 
affect the observed spectra In a manner that is differ¬ 
ent from the effect of changing the mass accretion rate, 
which can, therefore, be measured independently of these 
parameters from the normalization of the flux density. 

Since the ISCO depends on the spin and the deviation 
parameters £3, ai3, and 022 , the question of whether the 
spin and the deviation parameters can be measured in¬ 
dependently has to be analyzed more carefully. Since the 
emitted flux of a Kerr-like black hole only differs signif¬ 
icantly from the emitted flux of a Kerr black hole with 
the same ISCO radius if the latter black hole is spinning 
rapidly as can be seen from Fig. [2J one would expect 


that the observed flux density is very similar for Kerr 
and Kerr-llke black holes with the same ISCO as long 
as the Kerr black hole has a small to intermediate spin. 
This Is Indeed the case and the corresponding observed 
spectra are practically Indistinguishable. 

For high values of the spin, however, the observed spec¬ 
tra differ significantly. In Fig. [5l I Illustrate this for par¬ 
ticular choices of the deviation parameters 013 and 022 . 
Fig. [5] shows the observed number flux density for Kerr 
black holes with spins a = 0.9M and a = 0.95M and for 
Kerr-like black holes with smaller spins but a nonzero 
values of the parameters ai 3 or 022 chosen such that 
the ISCO coincides with the ISCO of the respective Kerr 
black hole. The observed spectra differ strongly at the 
high-energy end. Nonetheless, spectra of Kerr-like black 
holes with nonzero values of the deviation parameters are 
still very similar to the spectra of Kerr black holes with a 
different spin and, therefore, a different ISCO. However, 
if the spin of the Kerr black hole Is near maximal, spectra 
of Kerr-like black holes can extend to such high energies 
that they cannot originate from a maximally spinning 
Kerr black hole. 

These results imply that the deviation parameters cris 
and a 22 can be measured independently of the spin If 
the radius of the ISCO is ~ M and the magnitude of the 
deviation is of order unity. This also means that the lo¬ 
cation of the ISCO can be measured robustly as long as 
the ISCO radius is not too close to the black hole corre¬ 
sponding to Kerr black holes with small to Intermediate 
values of the spin. Since the limb darkening does not af¬ 
fect the high-energy end of the spectrum, nonzero values 
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FIG. 6 . Direct images of the inner parts of geometrically thin accretion disks around (left) a Kerr black hole with spin a = 0.95M 
and (right) a Kerr-like black hole with the same spin and a value of the deviation parameter ais = — 1 for fixed values of the 
inclination i = 80°, mass M = IOMq, distance D — IQ kpc, mass accretion rate M — 10^® g s“^, and spectral hardening factor 
/col = 1.7. Both panels show the observed number flux density at 1 keV in units of keV“^ cm“^ s“^. The highest emission 
originates from a strongly localized region near the ISCO on the side of the black hole that is approaching the observer where 
Doppler boosting and beaming are particularly high. This region shifts toward the black hole and emits a significantly higher 
flux for negative values of the parameter ais. 


of the parameter £3 can always be closely mimicked by a 
Kerr black hole with a different spin. 

In Fig. [ 6 l I illustrate the appearance of accretion disks 
for black holes with mass M = IOM 0 , spin a = 0.95M, 
and different values of the deviation parameter ai 3 as 
“seen” by a distant observer at an energy of 1 keV. The 
other parameters are held fixed with values of the incli¬ 
nation i = 80°, distance D = 1{) kpc, mass accretion rate 
M = 10^® g s“^, and spectral hardening factor fcoi = 1.7. 
Here I focus only on the innermost part of the disk with a 
radial extent stretching from the ISCO to r = 15M. This 
part of the disk experiences the most dramatic effects if 
the deviation parameters are nonzero. The observed flux 
density is strongly peaked on the side of the disk that ap¬ 
proaches the observer due to enhanced Doppler boosting 
and beaming. Decreasing the value of the parameter ai 3 
increases the peak flux drastically. Likewise, the ISCO 
radius decreases and the location of the peak flux shifts 
to a smaller radius. 

Because of their minute angular sizes, however, such 
images of accretions disks are unobservable for stellar- 
mass black holes. While the microquasar GRS 1915-1-105 
has a huge accretion disk with an estimated radius of 
^ 10^® cm that would be large enough to be imaged 
with a very-long baseline interferometer (see Ref. [ 63 ), 
its innermost region will not be resolved. In the case of 
supermassive black holes such as the one at the Galactic 
center, the observed spectra are much better described by 
radiatively inefficient accretion flow models (e.g., (H). 

V. RELATIVISTICALLY BROADENED IRON 
LINES 

In this section, I simulate relativistically broadened 
iron lines. As in Refs. [ 2 ^ [s^, I assume that the disk 
plasma moves on circular equatorial orbits at the local 
Keplerian velocity. I further assume that the emission 


from the disk is either isotropic or limb darkened as well 
as monochromatic with a rest frame energy Eq (e.g., 
Eq Ri 6.4 keV for the iron Ka line) and has an emissivity 
e(r) oc r““, where a is the emissivity index. 

The observed specific flux is then given by the expres¬ 
sion 

^ ^ j j - Eog{x',y')], (38) 

which I evaluate using the algorithm described in 
Refs. [ 2 ^, with the metric in Eq. m as the under¬ 
lying spacetime and the changes outlined in Section Hill 

The effect of the deviation parameters on observed 
iron line profiles is similar to their effects on the ther¬ 
mal disk spectra described in Sec. IIVBI In Fig. [71 I plot 
iron line profiles for black holes with spins a = QAM and 
a = 0.8M for several values of the deviation parameter 
022 . The disk inclination is i = 30°, the outer disk radius 
is r-Qut = lOOM, the emissivity index is a = 3, and the 
emission is isotropic. The energy E is displayed in units 
of the energy of emission Eq . For decreasing values of the 
parameter 022 , the fluxes of the “blue” and “red” peaks 
increase and the “red tail” of the line profile is shortened. 
The first effect is primarily caused by the orbital velocity 
of the accretion flow, while the second effect is deter¬ 
mined by the location of the ISCO and the photons that 
are emitted near the ISCO, which experience a strong 
gravitational redshift. 

The effect of the parameter ai 3 is very similar with 
the difference that the corresponding modification of the 
line profile as mentioned above occurs for increasing val¬ 
ues of the parameter ai 3 instead of decreasing values. 
Similar to the thermal disk spectra, the iron line profiles 
depend only marginally on the parameters £3 and a 52 if 
the emission is isotropic. 

At higher disk inclinations and for nonzero values of 
the parameters q;i 3 or a 22 , the blue peak is slightly al¬ 
tered, while the red peak is affected only marginally. Ref- 











11 




E/E„ 


FIG. 7. Iron line profiles for black holes with spins (top) 
a = 0.4M and (bottom) a = 0.8M with an outer disk radius 
Tout = lOOM, a disk inclination i = 30°, and an emissivity 
index a = 3 for several values of the deviation parameter 022 • 
The energy E is measured in units of the energy of emission 
Eq. The line profiles are altered primarily at high energies 
and in their extent toward low energies. 


erence [ 2 ^ observed a similar effect in the case of iron line 
profiles described by their Kerr-like metric which depends 
on only one deviation parameter [H and identified the 
flux ratio of the two peaks as a potential observable of 
a violation of the no-hair theorem. In the metric I use 
in this paper, this effect still prevails. However, for line 
profiles simulated in both spacetimes (as well as the Kerr 
metric) the red peak is often submerged into the line pro¬ 
file and may be difficult to identify in practice. Similar 
changes of the peak flux can also be achieved by different 
values of the emissivity index and the outer disk radius 
as is well known for Kerr black holes. These parameters 



E/Eo 


FIG. 8 . Iron line profiles for black holes with spin a = 0.8M, 
an outer disk radius rout = lOOM, a disk inclination i = 60°, 
and an emissivity index a = 3 for several values of the param¬ 
eter £ 3 . In all cases, the disk emission is limb darkened and 
the line energy E is measured in units of the emission energy 
Eq. For increasing values of the parameter £ 3 , the peak flux 
increases. 


have to be determined from a spectral fit of the entire 
line profile. 

If the emission is limb darkened, changing the parame¬ 
ters ai 3 and 022 has a similar effect on the line profiles as 
for Kerr black holes and the line profiles appear slightly 
narrower with an increased peak flux and, at lower incli¬ 
nations, a decreased flux at lower energies. For nonzero 
values of the parameters £3 and q; 52 , the line profiles re¬ 
main nearly unaffected at lower inclinations. At higher 
inclinations, however, the peak flux is modified and in¬ 
creases for positive values of the parameters £3 and 052 
and decreases for negative values of the parameters £3 and 
0 ^ 52 , which I illustrate in Fig. El Note that this modifica¬ 
tion is slightly different from the one caused by nonzero 
values of the deviation parameters Q!i 3 and a 22 for either 
isotropic or limb darkened emission. 

As in the case of Kerr black holes, for iron lines that 
originate from the Kerr-like compact objects that are de¬ 
scribed by the metric proposed by Ref. [13, Ref. HI] 
found that at least for small to intermediate disk incli¬ 
nations the inclination angle can be robustly determined 
from the location of the “blue edge” of the line, which 
depends only very little on the other system parameters 
including the deviation from the Kerr metric. For iron 
line profiles simulated in a background described by the 
metric in Eq. o, this is likewise the case. 

References [H, [53 also found that for arbitrary spins 
over the entire spin range iron lines of Kerr-like compact 
objects described by either their metric or the quasi-Kerr 
metric are strongly correlated with the iron lines of Kerr 
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FIG. 9. Iron line profiles for different values of the spin and the deviation parameters aia (top) and a 22 (bottom) such that for 
both sets of parameters in each panel the ISCO coincides. The other system parameters are held fixed with values rout = lOOM, 
i = 30°, and a = 3. For Kerr black holes with low to intermediate values of the spin, the line profiles are very similar to the 
profiles for Kerr-like black holes which have the same ISCO. For Kerr black holes with high spins, however, the line profiles are 
different from the profiles which correspond to the same ISCO radius with nonzero values of the deviation parameters aia and 
0:22 if the deviation is snfficiently large. 


black holes if the spin and deviation parameter are cho¬ 
sen such that the ISCO (or, more generally, the inner 
edge of the accretion disk [^) in both cases coincides. 
In the new Kerr-like metric, for nonzero values of the 
parameters 0:13 and 022 , such profiles are likewise practi¬ 
cally indistinguishable if the Kerr black hole is spinning 
slowly or moderately. For Kerr black holes with high 
spins, however, this strong correlation between the spin 
and the deviation parameters ai 3 and 022 does not per¬ 
sist and the line profiles can be significantly different from 
each other if the deviation is large enough. I illustrate 
this in Fig. |9l In practice, the emissivity index is a func¬ 
tion of the radius and can reach much larger values (^ 7 
or higher) near the ISCO. In this case, iron line profiles 
differ even at small values of the deviation parameters 
ai 3 and 022 , because the innermost part of the accretion 
disk where the deviations are most important due to the 
strong relativistic effects contribute significantly more to 
the total observed line flux. 

In their initial study of relativistically broadened iron 
lines in a non-Kerr background, Ref. [13] suggested 
that this strong correlation should indeed be broken for 
rapidly spinning black holes, because the accretion disk 
in this case extends almost all the way down to the event 
horizon, where the effect of a deviation from the Kerr 
metric is most apparent. While this did not occur in the 
metric of Ref. [l^ , the suggestion by Ref. [13| is indeed 
verified for deviations described by the parameters 0:13 
and Q !22 in the metric given by Eq. ([3]). Although it is 
difficult to compare these two metrics directly, because 
they are of two distinct forms, the reduced correlation 
seems to be caused by the different dependence of the 
respective metric elements on the deviation parameters. 
The two deviation parameters ai 3 and a 22 in the metric 
used in this paper affect only its (t,t), and {4>,4>) 


components. In contrast, the metric of Ref. [T^ depends 
on only one deviation parameter that affects all nonzero 
metric components apart from the {0^9) component. 

In addition, contrary to the case of the thermal disk 
spectra analyzed in Section fiVBi iron lines in this part of 
the parameter space cannot be mimicked by the profiles 
of black holes with a different spin, because this would 
manifest in an altered tail of the line at low energies due 
to the effect of the gravitational redshift. This illustrates 
an important difference between both methods: At high 
spins, thermal disk spectra already differ when the de¬ 
viation parameters ai 3 or 022 are small, while iron line 
profiles differ only if the deviations or the emissivity in¬ 
dex are sufficiently large. On the other hand, for sets 
of (high) spin values and deviation parameters ai 3 and 
a 22 with different ISCO radii, thermal spectra can still 
be very similar, while iron fines cannot be confused, at 
least in principle. 


VI. QUASIPERIODIC VARIABILITY 


In addition to the thermal accretion disk spectra and 
relativistically broadened iron lines, quasiperiodic vari¬ 
ability constitutes a potential third approach to measure 
the spins of black holes and even of potential deviations 
from the Kerr metric. In the following, I investigate the 
dependence of QPO frequencies on the spin and the de¬ 
viation parameters in two particular QPO models, the 
diskoseismology mod el ( see Ref. ^ for a review) and 
the resonance model [^, j^. Further viable QPO mod¬ 
els exist (e.g., El) , but I do not consider them here, 
because these do not predict fixed values of the spin (or 
the deviation parameters) at a given QPO frequency. 

In the diskoseismology model, QPOs can arise as so- 
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FIG. 10. Contours of constant g-mode frequency for a IOMq black hole as a function of the spin and the deviation parameters 
ai 3 (left) and 022 (right). In the top panel, frequency contours are shown as solid blue curves with frequencies (top to bottom) 
50 Hz, 100 Hz, 200 Hz, 300 Hz, 400 Hz, 500 Hz, 600 Hz, and 700 Hz. For reference, contours of constant ISCO radius are 
shown as dotted magenta lines with radii (top to bottom) 8M, 7M, ..., 2M. In the gray dotted region, the radial epicyclic 
frequency does not vanish at the ISCO and g-modes may not exist. In the bottom panel, frequency contours are shown as solid 
blue curves with frequencies (left to right) 75 Hz, 100 Hz, 200 Hz, 300 Hz, 400 Hz, 500 Hz, and 600 Hz. Contours of constant 
ISCO radius are shown as dotted magenta lines with radii (left to right) 6M, 5M, ..., 2M. In both cases, the contours of 
constant g-mode frequency are mostly aligned with the contours of constant ISCO radius except for larger deviations from the 
Kerr metric if the ISCO lies at a radius smaller than roughly 3M. The black shaded region marks the excluded part of the 
parameter space. 


called gravity modes (g-modes ( 3 ) and corrugation 
modes (c-modes 0 )- Here I focus on the fundamental 
g- and c-modes. In the resonance model, QPOs can be 
identified as resonances between the dynamical frequen¬ 
cies which preferentially occur at small-integer frequency 
ratios. As an example, I focus on a 3 : 1 resonance 
between the Keplerian and radial epicyclic frequencies; 
other resonances should have a very similar dependence 
of the deviation parameters. 

Expressions for the Keplerian as well as the radial 
and vertical epicyclic frequencies for particles on circular 
equatorial orbits in the metric given by Eq. can be 
found in Ref. [d^. The fundamental g-mode occurs at 
the radius where the radial epicyclic frequency reaches 
its maximum and the fundamental c-mode corresponds 
to the Lense-Thirring frequency (i.e., the difference be¬ 
tween the Keplerian and vertical epicyclic frequencies) at 
the ISCO [3)0- A parametric resonance between the 
dynamical frequencies usually occurs at a different ra¬ 
dius, which is determined directly by the frequency ratio 
of the resonance [13, |48| . 

In Eigs. [Mia I plot contours of the fundamental g- 
and c-modes and of the Keplerian frequency in the 3 : 
1 resonance, respectively. Eor reference, I likewise plot 
contours of constant ISCO radius. The contours of the g- 
mode and resonance frequencies are mostly aligned with 
the ISCO contours except for very high frequencies. The 


contours of the c-modes are less aligned with the ISCO 
contours especially for very low frequencies in the case of 
deviations described by the parameter ai^. 

The dependence of these QPO modes on the spin and 
the deviation parameters is similar to the one in the 
quasi-Kerr metric and the metric of Ref. dl as discussed 
in Refs. [H, [s^, [Ml. However, since at the lowest order in 
the deviation functions the metric in Eq. contains four 
deviation parameters instead of only one as in the other 
two metrics, the departure of the frequency contours from 
the contours of constant ISCO radius are generally larger. 

Since the dynamical frequencies and the location of the 
ISCO depend only marginally on the deviation parameter 
£3 and not at all on the deviation parameter 052 (apart 
from the radial epicyclic frequency, which is affected only 
slightly), their effect on the g-, c-, and resonance modes 
are small. 


VII. CONCLUSIONS 

In this paper, I analyzed three potential high-energy 
probes of deviations from the Kerr metric parametrized 
by a new Kerr-like metric : thermal continuum spec¬ 
tra, relativistically broadened iron lines, and quasiperi- 
odic variability. The Kerr-like metric describes a black 
hole and depends on four deviation parameters which 
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FIG. 11. Contours of constant c-mode frequency for a IOMq black hole as a function of the spin and the deviation parameters 
ai 3 (left) and 022 (right). In the top panel, frequency contours are shown as solid blue curves with frequencies (left to right) 
5 Hz, 10 Hz, 25 Hz, 50 Hz, 100 Hz, 200 Hz, and 400 Hz. For reference, contonrs of constant ISCO radins are shown as dotted 
magenta lines with radii (top to bottom) 8M, 7M, ..2M. The c-mode contours are only roughly aligned with the ISCO 
contours for large values of the spin and for deviations |ai 3 | < 3. In the gray dotted region, the radial epicyclic frequency does 
not vanish at the ISCO and c-modes may not exist. In the bottom panel, frequency contours are shown as solid blue curves 
with frequencies (left to right) 5 Hz, 10 Hz, 50 Hz, 100 Hz, 500 Hz, 1000 Hz, 2000 Hz, and 4000 Hz. Contours of constant 
ISCO radius are shown as dotted magenta lines with radii (left to right) 6M, 5M, ..., 2M. The contours of constant c-mode 
frequency are mostly aligned with the contours of constant ISCO radius except for larger deviations from the Kerr metric if 
the ISCO lies at a radius smaller than roughly 3M. The black shaded region marks the excluded part of the parameter space. 


affect different components of the metric. This work 
extends previous analyses that focused on other Kerr- 
like metrics 0 ,[I 1 which depend on only one deviation 
parameter and generally harbor naked singularities (see 
Ref. H). 

I showed that all three observables depend significantly 
on two deviation parameters (013 and 022 ) which mod¬ 
ify the and (^,0) elements. I also showed 

that in the case of limb darkenend emission a third devi¬ 
ation parameter, e^, which alters all nonvanishing metric 
elements affects both the thermal spectra and iron line 
profiles, while a fourth deviation parameter, 052 , which 
modifies only the (r, r) component affects only the iron 
line profiles at higher inclinations. The latter two param¬ 
eters have a relatively small effect on QPO frequencies 
and only lead to potentially detectable frequency shifts 
if their magnitudes are very large. 

Since the ISCO of the spacetime depends on the spin 
and the parameters 013 , 022 , and € 3 , it is important to 
characterize the non-Kerr effects on the three observ¬ 
ables for those parameter combinations that belong to 
the same ISCO radius. I showed that thermal spectra 
and iron line profiles are practically indistinguishable if 
the ISCO lies at a sufficiently large radius corresponding 
to small to intermediate spin values of Kerr black holes. 
For Kerr black holes with high spins, however, both the 
thermal spectra and iron line profiles differ for sets of 


parameters with the same ISCO radius. This difference 
manifests already for small values of these deviation pa¬ 
rameters in the case of the thermal spectra and for larger 
values of these parameters in the case of the iron lines un¬ 
less the emissivity index is very high. Such a difference is 
qualitatively new and has not been seen in other Kerr-like 
metrics (cf. e.g., HI). 

I also showed that in the diskoseismology and reso¬ 
nance models, QPO frequencies likewise depend strongly 
on the deviation parameters 013 and 022 and only 
marginally on the other two. In both models, contours of 
constant frequencies are mostly aligned with the contours 
of constant ISCO radius with exceptions at high spins 
and for low-frequency c-modes. In Table I, I summarize 
the effects of the deviation parameters on the various ob¬ 
servables. 

Already, a wealth of x-ray data has been collected 
which led to a number of spin measurements with the 
continuum fitting and iron line techniques (see, e.g.. 
Refs. [11,135). In order to constrain potential deviations 
from the Kerr metric, the ability to unambiguously iden¬ 
tify such deviations is essential. Unfortunately, all ob¬ 
servables discussed in this paper strongly correlate with 
the location of the ISCO at least in the case of black 
holes with small to intermediate spins. Therefore, these 
techniques will likely only be able to constrain combina¬ 
tions of values of the spin and the deviation parameters 
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FIG. 12. Contours of constant Keplerian frequency for a 3 : 1 resonance between the Keplerian and radial epicyclic frequencies 
for a IOM 0 black hole as a function of the spin and the deviation parameters 013 (left) and 022 (right). In the top panel, 
frequency contours are shown as solid blue curves with frequencies (top to bottom) 100 Hz, 150 Hz, 250 Hz, 450 Hz, 600 Hz, 
and 900 Hz. In the gray dotted region, the radial epicyclic frequency does not vanish at the ISCO and resonance modes may 
not exist. The additional blue line at the bottom forms the disjoint second branch of the 250 Hz contour which emerges at 
the boundary of the gray dotted region at a value of the spin a « —0.2M (not shown). For reference, contours of constant 
ISCO radius are shown as dotted magenta lines with radii (top to bottom) 8 M, 7M, .. ., 2M. In the bottom panel, frequency 
contours are shown as solid blue curves with frequencies (left to right) 200 Hz, 300 Hz, 500 Hz, 1000 Hz, and 1500 Hz. Contours 
of constant ISCO radius are shown as dotted magenta lines with radii (left to right) 6 M, 5M, ..., 2M. In both cases, the 
contours of constant y-mode frequency are mostly aligned with the contours of constant ISCO radius except for large values of 
the parameter 0 : 22 . The black shaded region marks the excluded part of the parameter space. 


Parameter ISCO Radius Thermal Spectrum Thermal Spectrum Iron Line Iron Line g-Mode c-Mode Resonance 


(isotropic) (limb darkened) (isotropic) (limb darkened) 


£3 

Moderate 

Weak 

Strong 

Weak 

Moderate® 

Weak 

Weak 

Weak 

ai3 

Strong 

Strong 

Strong 

Strong 

Strong 

Strong 

Strong 

Strong 

a22 

Strong 

Strong 

Strong 

Strong 

Strong 

Strong 

Strong 

Strong 

Oi52 

None 

Weak 

Weak 

Weak 

Moderate® 

Weak 

None 

Weak 


High disk inclinations only. 

TABLE I. Effects of the deviation parameters es, 0 : 13 , 022 , and 052 on the location of the ISCO as well as on thermal accretion 
disk spectra, relativistically broadened iron lines, and QPOs. The parameters 013 and Q 22 have a strong effect on all quantities. 
In the case of thermal disk spectra and iron lines, the effect of the parameters £3 and 052 depends on the emission type of the 
accretion disk (isotropic or limb darkened). The location of the ISCO radius is affected only slightly by the parameter £3 and 
is independent of the parameter 052 . These two parameters likewise have a marginal effect on the QPOs. 


which roughly correspond to the same ISCO radius un¬ 
less the black hole has a high spin. The consideration 
of additional deviation parameters, in the form of either 
higher order terms in the deviation functions given by 
Eqs. (H-dll) or generalized metrics which contain more 
than four such functions, would most likely only increase 
the level of degeneracy between spacetimes with similar 
ISCO radii but with different values of the spin and the 
deviation parameters. Nonetheless, an analysis of such 
observations in the Kerr-like background should be able 
to place first constraints on deviations from the Kerr met¬ 


ric. For some sources, two or even three of the above ob¬ 
servables have been detected and a combination of these 
measurements may further improve the constraints. 

Reference estimated the precision required of fu¬ 
ture x-ray missions such as Astro-H and ATHENA-I- to 
measure deviations from the Kerr metric parametrized by 
the metric of Ref. [l^ with iron lines. They found that 
black holes with higher spin values allow for tighter con¬ 
straints on the deviation parameter and that a precision 
of ~ 5% leads to constraints on the deviation parameter 
to less than order unity for spin values |a| ^ 0.5M. Since 
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FIG. 13. Rough constraints (light blue regions) on the spin a and deviation parameters ai 3 (left) and Q 22 (right) from the 
stellar-mass black hole Cygnus X-1. The light blue regions show the range of thermal spectra as a function of the spin and the 
respective deviation parameter for which the flux density at energies 0.7 keV < E < 4.0 keV lies between the flux densities of 
the spectra of Kerr black holes with spins a = 0.983M and a = M corresponding to the 3a limits of the recent spin measurement 
by Gou et al. [^. The other system parameters are held fixed at their values given in Ref. [^. In both panels, the allowed 
range of spectra is a narrow region that does not lie along lines of constant ISCO radius (blue dotted lines), because at these 
high spin values the relativistic effects in the innermost region of the accretion disk affect the spectra significantly. In the left 
panel, the spin and the parameter 013 are constrained to the approximate intervals a > 0 . 68 M and — 2.1 < 013 < 0 . 2 , while, in 
the right panel, the constrained region is bound only by the lower limit 022 > — 0.1 and extends to small values of the spin and 
very large values of the deviation parameter 022 (shown up to a value 022 = 3). The black long-dashed line labeled “Kerr/GR” 
corresponds to a Kerr black hole. The black shaded region in the left panel marks the excluded part of the parameter space. 


the parameters ai 3 and 0:22 in the metric of Ref. 
affect the line profiles in roughly the same manner (by 
magnitude) as the deviation parameter in the metric of 
Ref. [l^, these requirements should be similar if the for¬ 
mer metric is used as the underlying spacetime. 

In the case of thermal disk spectra, the Galactic black 
hole Cygnus X-1 is a prime candidate for a constraint 
on potential deviations from the Kerr metric. Assum¬ 
ing that this object is a Kerr black hole, Gou et al. 
recently measured a spin value a > 0.983M at Str confi¬ 
dence. In order to find a very rough estimate of such a 
constraint, I determined graphically the range of spectra 
as a function of the spin and the parameters 013 and 022 
for which the flux density lies between the flux densities 
of the spectra of Kerr black holes with spins a = 0.983M 
and a = M at energies 0.7 keV < E < 4.0 keV. In all 
cases, I considered only one deviation parameter and set 
the other deviation parameters to zero. The remaining 
system parameters were held fixed at their values given in 
Ref. [53 (inclination i = 27.1°, mass M = 14.8 Mq, dis¬ 
tance U = 1.86 kpc, spectral hardening factor /coi = 1.6). 
Since the normalization of the spectra is arbitrary in this 
analysis, I simply set the mass accretion rate to a con¬ 
stant value M = 10^® g s“^. All spectra assume limb 
darkened disk emission. 

This procedure results in two fairly narrow ranges of 


allowable values in the ( 0 , 013 ) and («, 0 ^ 22 ) parameter 
spaces, respectively. I plot these regions in Fig. [T31 Note 
that these regions do not lie along lines of constant ISGO 
radius, because the relativistic effects in the innermost 
part of the accretion disk are important at such high spin 
values. I obtain the approximate constraints a > 0.68M 
and — 2.1 < 013 < 0.2 in the case of deviations described 
by the parameter 013 as well as 022 > — 0.1 in the case 
of deviations described by the parameter 022 . In the 
latter case, the allowable range of spectra extends to low 
values of the spin and very large values of the deviation 
parameter 022 (not shown). 

A measurement of the parameters €3 or 052 with either 
thermal spectra or iron line profiles will be much more 
complicated, because their overall effect on the spectra 
is small. Other observables such as x-ray polarization 
[ 2 ^ 1 ^ a s well as time-resolved iron line profiles and time 
lags |66j| may constitute additional probes of deviations 
from the Kerr metric. 


A different approach is the imaging of the shadows of 
the Galatic center black hole Sgr A* and of the supermas- 
sive black hole in M87 with the Event Horizon Telescope 
671. Reference [2l| showed that the shapes of the shad¬ 


ows of supermassive black holes described by the metric 
in Eq. m depend only on the spin and the deviation pa¬ 
rameters ai 3 and 022 - Consequently, using either x-ray 
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or (sub)mm observations, deviations described by the pa¬ 
rameters ai 3 and 022 will probably be easier to detect in 
practice, at least in the near future. 


where p® is the (contravariant) photon 4-momentum, p® 
the photon 3-momentum, p^ its component normal to 
the disk, and is the 4-velocity of the disk plasma. 
Since from Eq. (ESI 
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Appendix A: PHOTON EMISSION ANGLE 


Pe = 


VW) 






(A3) 


n-pl 



(A4) 


Here, I derive an expression for the emission angle C 
of a given photon momentum vector relative to the disk 
normal. My derivation follows closely Refs. [s^. [^. 

The emission angle of the photon momentum vec¬ 
tor is measured in a local inertial frame with a basis 
{ep Cf, Cg, e^}. In the equatorial plane, where the ac¬ 
cretion disk is located, the disk normal h is simply e^, 
which is related to the coordinate basis {et,er,eg,e^} of 
the metric by the equation [2l| 

h = e.~ = -^=eg . (Al) 

The emission angle is, then, given by the expression 


From Eq. (l32ll . the redshift is given by the expression 


9 = 


obs 

Pt 


I 

m“P® ’ 


(A5) 


Therefore, the emission angle can be expressed as in Eq. 
(IMl) . In the case of a Kerr black hole, this expression 
reduces to the equation 


cosC = 


re 


(A6) 


cosC = 


Ip® I 


n°‘p 

u^^p 


e 

— 

e ’ 
a 


(A2) as in Refs. jH, li^. 
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